Slowly rotating black holes in alternative theories of gravity 
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We present, in closed analytic form, a general stationary, slowly rotating black hole, which is 
solution to a large class of alternative theories of gravity in four dimensions. In these theories, 
the Einstein-Hilbert action is supplemented by all possible quadratic, algebraic curvature invariants 
coupled to a scalar field. The solution is found as a deformation of the Schwarzschild metric in 
General Relativity. We explicitly derive the changes to the orbital frequency at the innermost stable 
circular orbit and at the light ring in closed form. These results could be useful when comparing 
General Relativity against alternative theories by (say) measurements of X-ray emission in accretion 
disks, or by stellar motion around supermassive black holes. When gravitational-wave astronomy 
comes into force, strong constraints on the coupling parameters can in principle be made. 
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I. INTRODUCTION 

General Relativity (GR) is an elegant theory which 
agrees with all observations at Solar System scale and 
beyond [J, i]; however its nonlinear, strong-field struc- 
ture still remains elusive and difficult to test Q. This, 
together with some long-standing problems in Einstein 
theory (like the presence of singularities, difficulties in 
explaining the accelerated universe and galaxy rotation 
curves, etc), has motivated the study of viable alternative 
theories of gravity. These theories, also known as modi- 
fied theories of gravity, aim to reproduce GR in the weak- 
field regime, but they can differ substantially from it in 
the strong curvature regime, where nonlinear effects be- 
come dominant. In order to pass current experiments, al- 
ternative theories should have the same post-Newtonian 
expansion as GR, at least to lowest order. However, large 
deviations are possible in relativistic systems: black holes 
(BHs). neutron stars, and cosmological models. 

BHs are natural candidates to investigate strong curva- 
ture corrections to GR. In the next decade, gravitational- 
wave detectors Q and high-frequency very long baseline 
interferometry (VLBI) [5| may provide direct observa- 
tions of these objects and of their nonlinear structure, 
completing the wealth of information from current elec- 
tromagnetic observations Q. The geometric structure 
of BHs encodes information about the underlying theory 
of gravity. Within GR, no-hair theorems (see Ref. 
and references therein) guarantee that stationary BHs 
are described by the Kerr solution and this assumption 
enters most of the calculations, including gravitational- 
wave emission, gravitational lensing and properties of the 
accretion disks. However, when corrections to GR are 
considered, BHs can support non-trivial hairs and new 
classes of solutions may exist. Hence, it is important to 
derive deformations to the Kerr metric arising from 

alternative theories of gravity and to predict astrophysi- 
cal observables within a more general, bias-independent 
framework. 



Previous studies on BH solutions in alternative theo- 
ries of gravity suffer from two majors limitations. First, 
given the plethora of alternative theories that have been 
recently proposed, most of the approaches have focused 
on a case-by-case analysis (with the notable exceptions of 
Refs. 0, [i3). Secondly, motivated and well-behaved cor- 
rections to GR are usually involved, so that BHs must be 
constructed numerically. In particular, rotating solutions 
are extremely challenging to find in closed form and the 
Kerr metric is usually regarded as unique in this context. 
Thus, analytical solutions describing rotating BHs in a 
broad class of alternative theories, as the one we present 
here, are of utmost importance. 



In this work, generalizing previous studies on static 
BHs we derive the metric of slowly rotating BHs 
arising as solutions of a large class of alternative the- 
ories of gravity, in which the Einstein-Hilbert action 
is supplemented by all quadratic, algebraic curvature 
terms coupled to a scalar field. Rotating BH solutions 
are relevant for several reasons. Astrophysical BHs are 
likely to be (rapidly) spinning, due to accretion effects. 
Thus, any realistic computation (for example the proper- 
ties of accretion disks) must take rotation into account. 
Furthermore, the imprints of possible strong curvature 
corrections are expected to be stronger for those pro- 
cesses taking place close to near-extremal rotating BHs, 
for which the curvature is larger. For example, the 
Kretschmann invariant, K, = Rabcd.R"''"^'^ , on the equa- 
torial event horizon of a Kerr BH of mass M and angu- 
lar momentum J = aM in Boyer-Lindquist coordinates 

reads K. = 48 Af^ [M + VAf ^ - a^] where here and in 
the rest of the paper we use G = c = 1 units. For a 
Schwarzschild BH (a = 0), /CM* = 3/4. However, for 
extremal Kerr BHs (a = M) this scalar invariant is ^ 60 
times larger, JCM"^ = 48. 
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II. GRAVITY WITH QUADRATIC 
CURVATURE CORRECTIONS 

We consider a class of alternative theories of gravity 
in four dimensions obtained by including all quadratic, 
algebraic curvature invariants, generically coupled to a 
single scalar field [ll|. The action of this theory reads 

^ =17- f V~9d^x [R - 2\/a<l>S/''<l> ~V {<!>) + h {4>)R^ 

lOTT J '- 
+ f2{<l>)RabR'"' + fsWRabcdR'"'"' + f M Rabcd* R""''"] 
+ S'mat [7(0)5^1^, *mat] , (1) 

where, in the matter action S'matj we have generically 
included a non-minimal coupling, which naturally arises 
in some string theories defined in the Einstein frame. In 
the following, we neglect the scalar self-potential V{(j)). 
Its inclusion, along with theories in asymptotically non- 
flat spacetimcs, is a natural extension of the present work. 

When /i = /z = -4/i and /a = /i, the the- 

ory reduces to the bosonic sector of heterotic string the- 
ory and the quadratic corrections reduce to the Gauss- 
Bonnet invariant. In that case matter is non-minimally 
coupled to gravity, 7(^) = e*^. Static BH solutions 
in Gauss-Bonnet gravity were found analytically in the 
small coupling limit |l2l [l3| and numerically for general 
coupling [3] (see also Ref. Stationary BHs with 

Gauss-Bonnet corrections were considered numerically in 
Ref. {T3| for slow rotations, whereas their highly spin- 
ning counterpart was recently constructed in Ref. [16|. 
Furthermore, when /i = /2 = /a = and = uaS, 
the above theory reduces to Chern-Simons gravity 
and slowly rotating BHs in this theory where obtained 
in Refs. [13, [3 • The field equations arising from Eq. ^ 
are explicitly given in Ref. where analytical, static 
BH solutions were also obtained in the small coupling 
limit. Here we generalize previous studies by construct- 
ing slowly rotating BHs in the general theory ([1]). 

The theory ([1]) has to be considered as an effective ac- 
tion, obtained as a truncation from a more general the- 
ory. For example in the low-energy expansion of some 



string theories, the Gauss-Bonnet and Chern-Simons 
terms arise as second order corrections in curvature. The 
Einstein-Hilbert term is considered as the first order term 
in a (possibly infinite) series expansion containing all 
possible curvature corrections. In this sense, GR may 
be only accurate up to O {aEj^^ and second order cor- 
rections may be important when dealing with nonlinear, 
relativistic solutions. For the same reason, we work in a 
pcrturbative regime in which possible higher order terms 
in H]) can be safely neglected. We consider the weak-field 
expansion of the coupling functions 

/,((/)) =?7,+a,0 + O(02), ^ = 1,2,3,4 

where rji and on are dimensionful coupling constants. 
When the coupling functions are constant, i.e. ai = 0, 
the theories above are usually labeled "non-dynamical" 
and they admit all vacuum GR solutions As a result, 
for small scalar fields the background solutions do not 
depend on rji. Although non-dynamical theories would 
have a different linear response, for example a different 
gravitational- wave emission [l^, [llj , here we are inter- 
ested in modified background solutions and we then focus 
on dynamical couplings. Remarkably, in the small cou- 
pling limit, the dynamical theory only depends on four 
couplings. Ui, regardless the coupling functions fi{(t>). 



III. SLOWLY ROTATING BLACK HOLES 

We consider the following metric ansatz for the sta- 
tionary, slowly rotating limit, 

ds^ = -/(r, e)df + g{r, e)~^dr^ - 2w(r) sin^ Odtdip + 
+r^eir, 6)d6^ + sin^ 6^{r, d)d(p^ (2) 

together with the scalar field (p = (f>{r, 6). In Appendix El 
we solve the field equations [ll[ order by order in a pcr- 
turbative scheme for slow rotations and small couplings. 
Here, we simply report the final result. We obtain that 
the slowly rotating BH metric functions read 
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whereas the scalar field reads 
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where the novel terms are those proportional to aa^ and 
to a^as in Eq. ([5]) and Eq. d?]), respectively. Interest- 
ingly, these terms are the dominant corrections at large 
distances, because they scale with a lower power of r 
than those proportional to 014. As explained in the ap- 
pendix, the metric is found by requiring asymptotic flat- 
ness and regularity for r > 0. The curvature invariants 
are regular in the exterior spacetime. The BH mass and 
angular momentum can be read off from the l/r coeffi- 
cients in Eqs. ^ and ([5]). The angular momentum reads 
J — aM, whereas the physical mass of the BH is Ai = 
M{1 + 49a|/(320M4)) The above solution is accu- 
rate up to order 0{a'^/M^,a1/M'^,aa^/M^) in the met- 
ric and up to order ©(a^/M^, a^JM^, aaf/M^, a^Ui/M^) 
in the scalar field. At this order, the angular metric func- 
tions O and $ are simply given by the slowly rotating 
Kerr solution. For a = 0, the slowly rotating BH cor- 
rectly reduces to the static one setting, in the notation 
of Ref. [lil, = a^/(167r), /3 = l/(47r), k ^ l/(167r). 
Furthermore, for as = 0, it reduces to the slowly rotat- 
ing Chern-Simons BH Slowly rotating black holes 
in Einstein-Dilaton-Gauss-Bonnet gravity can obtained 
from our solution setting 04 = and their exact metric 
is presented here for the first time. We have compared 
the analytical results with numerical solutions [15| , find- 



ing very good agreement. Interestingly, this solution only 
depends on the couplings as and 04, since the terms 
proportional to ai and a2 do not contribute to this or- 
der. Moreover, the corrections to the scalar field aris- 
ing from as and 04 enter at different order in a: the 
Kretschmann correction only introduces even powers of 
a, while the Chern-Simons term only introduces odd pow- 
ers [ll|. Nevertheless, both corrections affect the gravit- 
omagnetic part of the metric, for example giving a mod- 
ified frame-dragging effect. Finally, the corrections pro- 
portional to as scale with a lower power of r than those 
proportional to a4. Hence, they are expected to be dom- 
inant at large distances. 



IV. GEODESIC STRUCTURE 

Many interesting and potentially observable effects 
around astrophysical BHs ultimately depend on how par- 
ticles move in the region few Schwarzschild radii away 
from the event horizon. For example, the inner prop- 
erties of the accretion disk are strongly affected by the 
location of the innermost stable circular orbit (ISCO) 
and, in turn, by the geodesic structure of the underly- 



ing spacetime. Most of the computations assume that 
the spacetime is described by a Kerr BH. However, de- 
formed solutions arising in alternative theories would also 
affect particle motion, with potentially observable conse- 
quences. In the modified theories considered here, test- 
particles follow spacetime geodesies. This follows from 
the conservation of the stress-energy tensor, V'T^i/ = 0, 
which is guaranteed by the diffeomorphism invariant ac- 
tion ([T]). In many situations the geodesic motion of mas- 
sive and massless particles is enough to fully describe 
many effects of astrophysical interest. 

We consider the following matter action for a point-like 
particle 



dt 



(8) 



where m is the mass of the particle and 7((/)) is a possi- 
ble coupling function between the matter and the scalar 
field. For low-energy modifications from heterotic string 
theory, ^ = e"^ . In the small field limit, we may write 



7(0) = 1 + 460 + 0(02) 



(9) 



where 6 = for minimal coupling and h = 1/4 in heterotic 
string theory. We focus on equatorial motion {9 = 7r/2, 
^ = 0). The radial geodesic motion on the equatorial 
plane can be derived from the equation 



V{r) 



g fhE'^ - fL^ + 2jEL 
7^ I f + fh 



-Sj , (10) 



where j = —uj{r), h = ^{r)r^, and (5 = 0, 1 for massless 
and massive particles, respectively. Here E and L are 
the energy per unit of mass and the angular momentum 
per unit of mass of the orbiting particle, respectively. 
For circular orbits at r = rc, the corresponding values of 
E and L can be found by imposing V{rc, Ec, Lc) = = 
V'{rc, Ec, Lc) and, for (5 = 1, the ISCO location is defined 
through ^^"(riscO: -E'c, Lc) = 0. Finally, the frequency at 
the ISCO reads 



ISCO 



''ISCO 



/ - jEc/Lc 
hEc/Lc + j 



(11) 



In line with our approximation scheme, we expand the 
geodesic quantities around their Schwarzschild value, i.e. 

X = + X(i)a + X(2)a2 + x(3)as + X^^^aj 

+A:(^)aas + X'-^^aal + AT^^'a^as + A(^'aa^(12) 
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where X schematically denotes Tc, Ec and Lc- In general, 
the coupling b introduces lower order contributions, like 
those proportional to a^. This is due to the lower order 
dependence of the scalar field in Eq. ([7]). For the same 
reason, such corrections do not arise for terms propor- 
tional to a4, since the odd-parity correction to the scalar 
field vanishes on the equatorial plane. Substituting the 
expansion (|12p and solving order by order, we obtain the 
following ISCO location and the frequency at the ISCO, 
normalized by the physical mass M., 
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V 2099520 V6 755827201/; 

where we have kept only dominant terms in b and we 
are considering corotating orbits only. Counter-rotating 
orbits can be simply obtained by inverting the sign of 
a. The behavior of the ISCO frequency depends on sev- 
eral couplings. For 6 = 0, the dominant correction is 
0{a\) and contribute to increase the frequency. The 
first corrections proportional to the BH spin are 0{aa') 
and 0(aa\) and they contribute to lower the frequency. 
However, when a non-minimal coupling is turned on, its 
effect is dominant jisj . The ISCO frequency gets negative 
0{baz) corrections. Since this is the dominant effect, a 
decreasing of the ISCO frequency could be seen as a gen- 
eral signature of non-minimal couplings, regardless the 
relative strength of a, as and 04. 

The same procedure can be applied to null geodesies, 
which are the trajectories of massless particles. In this 
case, it is easy to show that the result does not depend 
on the coupling 7. We get 



DISCUSSION AND CONCLUSION 



We have found slowly rotating BHs, solutions of a class 
of alternative theories as general as the action ([T]). This 
theory supplements GR by all quadratic, algebraic cur- 
vature terms coupled to a scalar field. Our solution is 
presented in closed form up to some order in the angu- 
lar momentum and in the coupling parameters. To the 
same order, we discussed the most relevant properties 
of the equatorial geodesic motion, giving the ISCO and 
light-ring frequencies. 

With the analytical solution at hand, several exten- 
sions of the present work are possible. The properties 
of the (modified) accretion disk can be used to constrain 
the parameters of the theory Q. Furthermore, the study 
of the geodesic structure can be generalized to include 
non-equatorial orbits and an analysis similar to Ref. (23j 
can be performed. Another interesting issue is the lin- 
ear response of the slowly rotating BH. Strong curvature 
corrections to GR affect the linear stability analysis (T5| 
and the gravitational- wave emission [24 1. 

In addition, several extensions of the present solution 
are conceivable. First of all, going further in the ap- 
proximation scheme, up to order a'af, corrections to the 
event horizon location and to the ergoregion would ap- 
pear. This can have a profound impact on the stability 
of these solutions. Furthermore, highly spinning BHs are 
phenomenologically more relevant and larger deviations 
from the Kerr metric may be expected. However, they 
have to be constructed numerically [iBl on the basis of 
a case-by-case analysis. In this case, our analytical so- 
lution can be useful; for example it can be used as an 
initial profile to start numerical relaxation methods, or 
to check numerical solutions. 
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where finuU = inuii/^'nuU is the light-ring frequency, and 
it is related to the real part of the ringdown frequency of 
the BH in the eikonal limit [22j. The dominant correc- 
tion is 0{a\) and contributes to increase the frequency, 
whereas the 0(aa§) and 0{aa\) corrections have an op- 
posite relative sign. 



We report here that the slowly rotating metric we 
found can be mapped into the bumpy BH formalism 
along the same line discussed in Ref. although the 
mapping is non-trivial. On the other hand, this solu- 
tion does not bel ong to the class of deformed Kerr BHs 
proposed in Ref. [lO|. 
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Appendix A: Slowly rotating approximation 

The modified field equations are obtained by varying 
the action ^ with respect to the metric and to the scalar 
field. Varying the action ([T]) with respect to the metric, 
neglecting the Smat term and the potential V{(j)), we find 

Gab + aiUab + Ct2lab + OL^J ah + OL^^ab = T'ab' (^1) 

where Tfj, = \/a(f>^b(f> ~ ^SahVct^V^ and Hab, ^ab, Jab, 
JCab are explicitly given in Ref. Varying the ac- 

tion ^ with respect to the scalar field 0, we get 

-2n^ = aiR'+a2RabR''''+azRabcdR'''""'+aiRabcdR''''"' 

We shall neglect terms of order in the equation 
above. Since the Ricci scalar and the Ricci tensor are 
both zero in the background spacetime, the scalar field 
equation reduces to 

n<j, = -i(«3i?,fc,rfi?°''-'' + «4 *RabcdR^''''), (A2) 

where the tilde stands for background quantities. We 
note here that, when a = 0, we recover the scalar 
field for spherically symmetric Gauss-Bonnet BHs jl^ . 
since at this order the Gauss-Bonnet term is just the 
Kretschmann invariant and there is no correction from 
the Chern-Simons term [l^. On the other hand, for 
Q!3 = we recover the scalar field for slowly rotating 
Chern-Simons BHs. Also, there is no correction of order 
a^a, since the Kretschmann invariant has only correc- 
tions in even powers of a. Therefore, we can write 

cf^^C^GBfiS ^^^^2^c^^^Q-^^ (A3) 



where (fpBfiS ^^^^ scalar field for spherically symmetric 
Gauss-Bonnet BHs plus the correction of slowly rotating 
Chern-Simons BHs, both assuming small coupling con- 
stants. Substituting Eq. (|X3|) in Eq. we find that 
the only solution for (j)'^ which is regular at the horizon 
and goes to zero in the limit ^ ^ 1, is given by the 
corresponding term in Eq. ([T]) . 

Considering corrections up to af, the modified Ein- 
stein's equations read 

Gab + 8a3^ahcdV^V'^(/) + 8a4*i?(a%)''VrfV,0 

= i (2Va0Vfc(/) - gabVc4>'^'<t) , (A4) 

in which the scalar field (f) is given by Eq. ([7]). We note 
here that the lowest dynamical corrections to the metric 
are given by second order terms in 03 and a^. Therefore, 
we can write 



GB,CS ,2c /Ar-N 

9ab ^ 9ab + "a^fab : (A5) 



where g^j, ' is the metrics for the spherically sym- 
metric Gauss-Bonnet BH plus the correction for slowly 
rotating Chern-Simons BH, both assuming small cou- 
pling constants. In the slowly rotating regime, the only 
non- vanishing term in g^j^ is gf^ With the ansatz 

g^^ ~ ~H{r) sin^ 6, we find that the only solution for 
H(r) that goes to zero in the regime ^ ^ 1, is given by 
the corresponding term in Eq. ([S])- 
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